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Abstract: 

We argue that a large fraction of the proton spin comes from a contribution due to a 
nonperturbative intrinsic charm component of the proton. An account of this contribution 
removes an apparent contradiction between the data and exact Kilim-Zakharov low energy 
theorem. On the other hand, we show that a large intrinsic charm spin component of the 
proton is implied by recent CLEO data on B —> r( decays. We argue that the proton 
spin and B —> rf data are manifestations of the same physics with full agreement between 
two different estimates of the intrinsic charm component of the proton spin. Our results 
suggest an explanation of the polarized DIS data and have profound consequences for 
future experiments on charm production off a polarized proton. A microscopic origin of 
the effect is related to strong well-localized gluon fluctuations (instantons) where the mass 
of the charmed quark is not a suppression factor at all, as was naively expected earlier. 



1 Introduction 


The famous proton spin problem is now almost ten years old. A wide interest in the 
structure of the proton spin was triggered by the 1987 EMC experiment |T| which has 
measured the first moment of the proton spin structure function Tf = / dxgl(x) = 0.126± 
0.018. This result implies that the polarized deep inelastic scattering (DIS) data cannot be 
explained simply by the valence quark spin components |2|, but, in addition, a substantial 
sea-quark polarization is needed. Such a polarized quark sea seems very miraculous from 
the point of view of both the relativistic quark model (where valence quarks account for 
3/4 of the proton spin) and the parton model (where the hclicity conservation prohibits 
a perturbative sea polarization by hard gluons). This problem, sometimes called “the 
proton spin crisis”, has stimulated a great deal of interest over the last decade (see e.g. 
P| for a review and representative list of references). 

Nowadays, it is commonly accepted that these are gluons that play an important 
role in the spin structure of the nucleon. Depending on a factorization scheme, a gluon 
contribution to T^ is either given by an independent anomaly term in the chiral-invariant 
scheme f|, or included in spin-dependent quark distributions in the Operator Product 
Expansion (OPE) based gauge-invariant scheme advocated by Jaffe and Manohar j5J. In 
the Jaffe and Manohar approach, which we will follow, the data are usually interpreted 
as the manifestation of a negatively polarized s-quark component in the proton. An 
equivalence of these two descriptions for the first moment T/ was established by Bodwin 
and Qiu |6|]. While the former presumably makes more contact with one’s intuition and 
the parton model, the latter is based entirely on QCD and turns out more convenient 
technically (see the Cheng’s paper in || and below). 

Irrespective of the question of interpretation of rj*, the gluon polarization A G(x) in 
the proton is a very interesting object on its own. This quantity is currently attracting a 
great deal of attention of both theorists and experimentalists. Among experimental tests 
on A G(x) a hadronic heavy-quark production (and, in particular, the charm production) 
is believed to be the best one. The reason for this belief is that within the standard 
perturbative photon-gluon fusion prescription the charm production amplitude is very 
sensitive to the gluon distribution in the nucleon. Among future experiments, intended 
to probe the gluon polarization through this scenario, one should mention the RHIC and 
HERA-A projects on a charm production in polarized pp collisions, and the COMPASS 
experiment at LHC on an open charm production in polarized DIS (see the review by 
Ramsey in || for references and more detail). Existing theoretical expectations for these 
experiments are all based on the perturbative photon-gluon fusion scenario. 

However, while there is little doubt that this mechanism is operative for unpolarized 
DIS, its relevance may (and must) be questioned in the case of polarized DIS. Indeed, 
in the latter case one deals with the axial channel instead of the vector one, as it was 
for usual DIS. On the other hand, it is very well known that physics of axial channels is 
drastically different from that of vector ones. Unlike the latter, axial and pseudoscalar 
channels are strongly influenced by nonperturbative effects. The best known 
example of this difference is provided by a nonperturbative breakdown of the Zweig rule 
in the pseudoscalar nonet. Thus, the perturbative charm production mechanism, valid 
for usual DIS, certainly cannot be taken for granted in the polarized case. In particular, 
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strong nonperturbative fluctuations in the axial channel could induce a sizeable “intrinsic 
charm” fraction of the proton spin, thus substantially changing theoretical predictions for 
charm production in polarized experiments. 

The purpose of this paper is to argue that the latter hypothetical situation is in fact 
precisely what happens in reality. We will present two independent estimates of the intrin¬ 
sic charm component of the proton spin, which agree with each other and both indicate 
that this quantity is quite large. Our main argument is based on recent experimental 
results, while a second estimate relies on a theoretical result which has been known for 
a long time, but unfortunately has not received much attention in the literature. As will 
be shown below, the observation of a large c-quark contribution to the proton spin turns 
out to be crucial for understanding the polarized DIS data. Our arguments go beyond 
perturbation theory. 

Our main estimate comes from recent data of B-physics (!), which at first sight might 
seem entirely uncorrelated with the proton spin problem. We mean the new CLEO results 
on B —> rj' decays 0. Recently, we have shown || that rj 1 production in B decays can only 
be explained by the presence of a large intrinsic charm component ( 0 | 07 ^ 750 )?/) of the ?/. 
Moreover, a theoretical prediction for this quantity j 8 | yields a good agreement with the 
data in both the exclusive [§] and inclusive 0 cases. On the other hand, it is fairly obvious 
that this property of the rj' implies a sizeable nonperturbatively generated intrinsic charm 
component of the proton spin (p\c^ fJ ,'y 5 c\p), as soon as there exists a Goldberger-Treiman 
type relation between these matrix elements. A relation with the R-physics provides us 
with an estimate of this quantity which, as will be shown below, agrees with an estimate 
based on a pure theoretical result which we are about to discuss. 

Our second estimate is based on a comparison of the DIS data with a low energy 
theorem due to Kuhn and Zakharov (KZ) H, which is an exact (in the chiral limit) 
formula for the proton matrix element of the topological density. This important result, 
obtained in 1990, sank into oblivion in subsequent publications. We understand that one 
of the reasons for this attitude was an apparent sign contradiction between the KZ formula 
and the data. As we will argue, this contradiction is actually removed by the intrinsic 
charm component of the proton spin in a very natural way and in full agreement with our 
previous estimate based on the analysis of B decays. Therefore, the two independent lines 
of reasoning lead to the same conclusion on a role of the polarized charm in the proton 
spin. 

We thus arrive at a coherent and attractive explanation of the polarized DIS phe¬ 
nomenology. In our opinion, the measured first moment Tf is determined (we assume the 
chiral limit) by four cornerstones: 

(1) spontaneous breaking of the SUl{3) x SUr(3) chiral invariance (see below), 

( 2 ) a resolution of the 17(1) problem, 

(3) the proton matrix element of the topological density, fixed by the KZ theorem, and 

(4) the intrinsic charm component of the proton spin, which is extracted from data on 
B —> 7 / decays. 

As we will show below, a consistent treatment of these four ingredients leaves very little 
(if any) of the mystery of T^, and has profound consequences for the future experiments. 
Our presentation is organized as follows. We start in Sect.2 with a re-interpretation of 
experimental facts on the polarized DIS. The only difference from the standard analysis is 


2 



that we assume a non-zero (but unspecified any further at this stage) value of the matrix 
element {p\cq^^c\p) which describes the intrinsic charm component of the proton spin. 
We use the data to select a 51/(3) flavor singlet contribution (which is a combination of 
the anomaly and intrinsic charm terms) to different flavor components of the proton spin. 
We thus interpret the data as a constraint on this 51/(3) singlet contribution. In Sect.3 
we present the KZ theorem and then use it to select from this combination a contribution 
of the intrinsic charm proper. In Sect.4 the same quantity is estimated independently 
from the data on B —> rj' decays. In Sect.5 we briefly discuss consequences of our results 
for a charm production in polarized experiments. Our conclusions are presented in final 
Sect. 6. 


2 SU(3) singlet proton spin from the data 

We start with the OPE analysis of the first moment of the polarized g\ structure function. 
The only deviation from the standard treatment discussed in detail in || is that we keep 
the charm component of the electromagnetic current of virtual photon. Retaining this 
term amounts to keeping trace of an intrinsic charm component of the proton spin. In the 
standard analysis this component is discarded on the grounds that this contribution is 
expected to appear only at higher orders in a s from loop effects. As will be shown below, 
this term is actually 0(a°), though it is suppressed by 1 /m 2 in agreement with general 
arguments by Jaffe and Manohar j5|. Therefore, we insist on keeping it in the analysis. 
On the other hand, at EMC energies Q 2 ~ 10 GeV 2 , the mass of a struck c-quark can be 
neglected, and thus the sole effect of taking the charm contribution into account in the 
OPE analysis^] is simply reduced to adding a corresponding charm contribution to light 
u, d, s quark spin components. In the standard nomenclature, we calculate the quantity 

r?« 2 ) = /' d xgi (x,Q 2 ) . (i) 

Jo 

Using the OPE for a T-product of two electromagnetic currents and selecting an antisym¬ 
metric in Lorentz indices contribution, one arrives at the result 

U« 2 ) = \ (5A«(Q 2 ) + lAd(Q 2 ) + 1 a S (Q 2 ) + jAcIQ 2 )) (1 - GA 2 + ...). (2) 

Here A q ( where q is any of the quark flavors) stands for the spin component of the proton 
due to the quark flavor q 

s^AqiQ 2 ) = (p\q^^ 5 q\p) , (3) 

(s M is the proton spin vector) and the operator in (|j) is normalized at the normalization 
point Q 2 . To select a SU( 3) singlet contribution to (|2]), it is convenient to use the 
combinations 

g\ = A u — Ad , g\ = A u + Ad — 2 As , (4) 

which have no anomalous dimensions and thus can be found from low energy neutron 
and hyperon beta decays data. Assuming the 577(3) flavor symmetry, the non-singlet 

Wsual unpolarized DIS on intrinsic charm was considered in ]lT|. 
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couplings are 


g\ = F + D , g\ = 3F — D . (5) 

For the values F = 0.459 ± 0.008 , D = 0.798 ± 0.008, this yields g\ = 0.579 ± 0.025 
[7|. Using relations (^), one can now express the answer for Tf in terms of the (known) 
constants g A ,g A and the unknown SU(3) singlet combination 

AE(Q 2 ) = Au(Q 2 ) + Ad(Q 2 ) + As(Q 2 ) (6) 

plus an intrinsic charm part (p|c 7 M 75 c|p). We obtain 

r; = C ns (Q 2 ) (Uj3 + + ic s (Q 2 ) (S(Q-) + 2Ac(Q 2 )) . (7) 

As we have said before, the only difference from the standard formulas in Eq.(j7j) is the 
presence of an intrinsic charm term Ac. Coefficients Cns, Cs account for perturbative 
QCD corrections Cns(s) = 1 — ot s /n + ... which can be found in [j|]. 

A global fit to all available data (including the data of SMC, E142 and E143) together 
with a treatment of higher order QCD corrections, but neglecting the charm contribution 
givesQ [0 


Aw = 0.83 ± 0.03 , Ad = -0.43 ± 0.03 , As = -0.10 ± 0.03 (8) 

with 

AS = 0.31 ±0.07 (9) 

at Q 2 = 10 GeV 2 . In our case a difference from this result comes due to retaining the 
intrinsic charm Ac term in Eq.((7). Therefore, we translate Eq.(|S]) in the constraint 

g° A = AE + 2Ac = 0.31 ± 0.07 . (10) 

(again Q 2 = 10 GeV 2 is implied). It is convenient at this stage to proceed to derivatives 
instead of axial currents which stand in (|10|) . Using the result of the OPE for the c-quark 
bilinear (see Eq.(|0) below) and working in the chiral SU(3) limit (which will be implied 
in what follows), we obtain for rif — 3 light flavors 

(p\ n f ^GG - 2 • j^—^g 3 GGG\p) = g° A ■ 2M p pi^ b p . (11) 

Here the first term is the standard contribution due to the light quarks. The second 
term is new, and plays a very important role in what follows. It originates in the c-quark 
contribution Ac, and technically comes from the derivative of the axial charmed current 

<9 m (c 7 m 75 c) = -^G“„G a ^ + 2m c ci'y 5 c (12) 

_ fia _ /~ia /ba _ ^ ^3 fcibc^a /~<b _i_ 

— Ur ^ 4vr 167r 2 m 2 ^ ^ U ' ’ ‘ ‘ 

2 These numbers were slightly updated by Ellis and Karliner in ||. To our accuracy this difference is 

unessential. 
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(see the appendix in || for a detailed derivation of this result). 

If we were to forget about the second charm-induced term in Eq.flTTD, we would con¬ 
clude that the quark contribution to the proton spin is given by the proton matrix element 
of the topological density (which has to be understood as the limit of a near-forward ma¬ 
trix element for vanishing momentum transfer )[]. The corresponding matrix element is 
fixed by a beautiful result due to Kuhn and Zakharov (KZ) [TCj], which provides us with 
an exact (in the chiral limit) answer for the matrix element of the topological density: 


ry _ 9r? f 

(p\n f ^-GG\p) = — ^2M p pijrj). (13) 

As will be shown below, this elegant formula is not sufficient to explain the polarized DIS 
data. The second charm-induced term in Eq.(|ll|) turns out to be absolutely crucial in this 
problem. We postpone a corresponding discussion to the next section, while here we want 
to go one step further and ask the question whether it is possible to fold flavor contribu¬ 
tions to the proton spin separately. The answer to this question is affirmative. To make 
contact with previous analyses which have not taken into account the charm contribution, 
we fold it convenient to split it equally between the three light flavor contributions, i.e. to 
consider “shifted” values A q' = Ag + 2/3Ac for q — u,d, s. This is in accordance with our 
definition ([H]) where Ac simply redefines the singlet contribution AS. Such a procedure 
amounts to adding an SU(3 ) singlet piece equally to each of the spin dependent light 
quark distributions, which was accounted previously by the anomaly term alone without 
a charmed contribution. We may then translate the result of Ellis-Karliner fit (|8]) in the 
values of A q' rather than A q as it was originally stated in Eq.@. On the other hand, 
taking the derivative for each of the three flavor q = u,d, s, we obtain 


(p\2m g qiy 5 q + Q\p) = A q' ■ 2M p piy 5 p , (14) 

where Q stands for the SU( 3) singlet part : 


Q 


T, Gd -\w^ GdG - 


(15) 


Here we would like to remind the reader that we work in the chiral limit m q = 0. However, 
it would be completely erroneous to put m q = 0 in matrix elements (|T4D from the very 


3 While, as will be explained below, this neglect is not legitimate, our point here is different. The fact 
that the light quark fraction of the proton spin A£ is equivalent to a particular value of the anomaly 
matrix element over the nucleon was strongly emphasized by Jaffe and Manohar Js|. This statement is a 
consequence of the exact equation of motion. A very similar situation holds in another anomalous case: 
the rj’ residue / v i measures a coupling of the r/ to quarks, ( 0 | W) = iVSfv'Qv On the other 

hand, the same parameter gives in the chiral limit a coupling of the 77' to gluons: (0|n/a s /(47r)GG|7/) = 
V 3 /„ '777^, as a consequence of exact equation of motion. As was stressed in j^], gluon operators cannot 
appear explicitly in the OPE (^]) since there is no gauge invariant gluon operator of spin 1 and twist 
2 . Thus, in the gauge invariant OPE approach to polarized DIS a gluon contribution is “hidden” (via 
the subtraction of regulator terms) in quark spin dependent distributions, which therefore do not have a 
parton interpretation. Later, Bodwin and Qiu have demonstrated ](J that any gauge invariant (i.e. Pauli- 
Villars) regularization of the photon-gluon scattering diagram automatically yields zero contribution of 
this diagram to the first moment in agreement with the general logic of OPE. Moreover, the approach 
based on the OPE is more convenient technically than the alternative method [Q, since it deals at every 
step with matrix elements of local gauge invariant operators, which are well defined objects. 
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beginning. Proceeding this way, we would miss the phenomenon of spontaneous breaking 
of the chiral SU(3)l x SU(3)r invariance. Goldstone bosons, resulting from this break¬ 
down, have masses m 2 ~ rn q , and therefore are the only surviving intermediate states 
in matrix elements {p\'m q qi'y^q\p) in the chiral limit m q —> 0, which is only taken after a 
saturation of these matrix elements by goldstones. Technically, it is more convenient not 
to do it explicitly, but rather use relations (d). Taking the derivatives, we obtain 

(p\2m u ui^ 5 u - 2m d di^d\p) = g :i A ■ 2M p pij 5 p , 

{p\2m u ui^ b u + 2m d dip b d - Am s si^s\p) = <f A ■ 2M p pi-/ 5 p , ( 16 ) 

(p\ 2 m u ui^ 5 U + 2m d di^d + 2m s si^s\p) = 0 . 


We would like to stress that the last of Eqs. fllGl) implies a resolution of the U(l) problem, 
i.e. the absence of a SU{3 ) flavor singlet meson with m 2 ~ m q . We are then able to find 
exact expressions for the SU{3) flavor dependent parts for each of the light flavor spin 
distributions separately: 


( p\2m u ui^u\p) 
{p\2m d di^ b d\p) 
(p\2m s sip 5 s\p) 


{^9a + ^ g \) • 2M v pi^p = (0.725) • 2M v pi^p , 

{~\dA + ^ 9a ) ■ 2M p p'hf 5 p = (-0.532) • 2M p pi- f5 p , (17) 

' 2 M pPHoP = (-0.193) • 2M p pi^ b p . 


In the numerical values in (|T7| ) we have implied that the numbers g\, g\ in the chiral limit 
are not very different from their phenomenological values in the real world (d) . We expect 
an accuracy of this approximation to be rather goodQ It is curious to note that, despite 
of their very transparent meaning and exactness in the chiral limit, we were unable to find 
relations ( |T7| ) in the literature. Meanwhile, they reveal something very interesting. Indeed, 
comparing these numbers to Ellis-Karlincr fit of experimental results (|8|) (remember that 
we have agreed to understand them as “shifted” values A q’ = A q + 2/3Ac), we see that, 
taken separately, each light flavor spin distribution is given by a large SU{3) dependent 
part plus a small SU{3 ) singlet piece whose numerical value is approximately the same 
and is very close to 0.1 for all the light flavors. We thus see that the mystery of 
T/ has to a large extent disappeared: the SU(3) flavor dependent parts give the main 
contributions to each of the light flavor spin distributions A q, but, on the other hand, these 
SU(3) variant parts cancel out in the total light flavor contribution AS, as a consequence 
of resolution of the U(l) problem (see the last of Eqs.(P^)). 

Therefore, if we managed to establish theoretically the number 


1 

2M p pij 5 p 




2 1 
3 167r 2 'm^ 


g 3 GGG\p) 


~ 0.1 , 


(18) 


(which is simply 1/3 g% see Eq.(|TT|)), the T^ problem would be resolved in QCD terms. 
Note that constraint (pA|), resulting from comparison of numbers ([T7| ) with Ellis-Karliner 

4 More accurately, it can be shown that a decrease of M p due to diminishing of the s-quark mass 
nearly compensates an increase of in the chiral limit. Thus, their product remains approximately 

the same as in the real world. 
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fit (||), agrees with Eq. (ED obtained without a specification of different flavor contribu¬ 


tions to the quantity £ + 2Ac. Thus, 577(3) singlet contribution (p~8[) is indeed flavor 
independent as should be expected. Furthermore, it consists of two parts. Had we knew 
the anomaly matrix element over the proton independently, we would then be able to 
extract the intrinsic charm component to the proton spin (the second term in Eq.(p^)) 
from experimental constraint (|l^). As will be shown in the next section, this knowledge 
is provided by the KZ theorem, and thus such a separation can indeed be made. 


3 Kuhn-Zakharov low energy theorem and intrinsic 
charm in the proton spin 

Our aim in this section is to separate the anomaly and intrinsic charm contributions in 
Eq. (|18]) . To this end, we are going to a use a beautiful result by Kuhn and Zakharov (KZ) 
m which provides us with an exact (in the chiral limit) answer for the proton matrix 


element of the topological density (the first term in (|l|)). KZ have shown that the 
famous dimensional transmutation phenomenon, which lies at heart of QCD low energy 
theorems, actually also fixes a value of the latter matrix element. In view of the fact that 
this important result is usually either ignored or doubted in the literature, we would like 
first to present KZ arguments at length and advocate their correctness. Next, we will use 
this result to obtain our first estimate of the intrinsic charm contribution to the proton 
spin. 

We start with recalling the well known fact that in a massless asymptotically free 
theory (such as QCD in the chiral limit) the appearance of a mass scale is related to 
the so-called dimensional transmutation phenomenon, according to which in this case the 
only mass parameter in the theory is 


m = A exp 


8vr 2 


(19) 


WW) 

where b = 11/3N C — 2/3 nj is the first coefficient of the Gell-Mann - Low beta function 
and A stands for a ultraviolet cut-off ( N c is a number of colors). By construction, m does 
not depend on the cut-off A, and all physical parameters (masses, vacuum condensates, 


etc.) in the chiral limit can only be proportional to a power of the mass parameter (pj9|) , 


as they cannot depend on the cut-off explicitly. Specifying on the proton mass M p ~ m, 
this means that under the variation of the cut-off 


a total variation of M p must vanish: 


A —>■ A(1 + e) 


eA d -^ = 0 

dA 


( 20 ) 


( 21 ) 


On the other hand, the total derivative in respect to A can be represented as the sum 
of a partial derivative and a term corresponding to a variation of the Lagrangian under 


cut-off shifts (j20|) : 


dM p dM p 

0 = eA ,. = eA^r .—b S\Ad p , 


dA 


da 


( 22 ) 
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where 


5 A M p pp = -(p\6 a L (A] \p) = -(psA 


OLA) 

dA 


I P) 


(23) 


and L^ is a regularized QCD Lagrangian. Its variation under scale transformations 
m is expressed by the conformal anomaly equation through the trace of the energy- 


momentum tensor 9 




f) r( A ) ho 

£d ^ - £ ^ Ga ^ • 

On the other hand, the partial derivative of M p is fixed by Eq. 

. dM v 

eA— = sM r . 


Combining Eqs. (|22|, ^3 j|24| , ^5|) , we obtain 


hfy 

(pI - b> = m pPp 


(24) 


(25) 


(26) 


Eq. (p6|) is well known and expresses the fact that in the chiral limit the proton mass is 
given by the conformal anomaly matrix element over the proton. As was pointed out in 
TOU , Eq. ([26|) can be obtained in a number of ways, thus providing an independent check 


of the above procedure. 

A particularly interesting observation made by KZ JT(J was that the very same line 
of reasoning could be applied in situations when only a fermionic cut-off (specifically the 
mass Mr of a Pauli-Villars fermion regulator) is varied in two possible ways: 


Mr —> Mr( 1 + £) 


(27) 


or 

Mr —■> M r ( 1 + ie ) , (28) 

where e is a small and real number. Consider first variation m- Again, the renor- 
malizability of the theory implies that a total variation of the proton mass under shifts 
(^7|) must vanish. On the other hand, a variation of the Lagrangian under a real shift 
of Mr can be readily found. Indeed, at the one-loop level contributions of fermion and 
boson regulators to the conformal anomaly are additive, and therefore the variation of 
the Lagrangian under transformations (0) is 


5m r L — £ ■ -rifOt. s G" , 


which together with Eq.(|2B|) yields 


, , dM p 2rif , r 
eMr ——— = — £ — -M p . 
BMr 3 b p 


(29) 


(30) 


Here we come to a most interesting (and sometimes wrongly criticized in the literature) 
part of the KZ proposal. Consider now the complex regulator mass variation (|28|) . This 





















variation leads to a y 5 mass term for the Pauli-Villars fermion R, as Eg. (|58|) is equivalent 
to the transformation 

MrRR —■> MrRR + £M R Ri r y 5 R . (31) 

Therefore, the variation of the Lagrangian in this case is proportional to the axial anomaly 


a* 


5iM R L — enf—GG . 


(32) 


Again, a shift in the physical proton mass due to the matrix element of anomaly (|32|) must 
be cancelled by a partial derivative in respect to variation (^8|) (see a comment below). 
On the other hand, the partial derivative is now determined by Eq.(j3(j) : 


is Mr 


dM ,, 


. 2 lit 


dM, 


p = -ie^±M v . 
r 3 b p 


(33) 


Therefore, matrix element of the chiral anomaly (^) is fixed 
independence of the physical mass of shifts (^) : 


T0j by the requirement of 


a ~ 2 

(p\^ gg \p) = , 


(34) 


that completes the proof of the KZ theorem. As was stressed in [TO] , the above derivation 


implies (again) a resolution of the U(l) problem since otherwise the matrix element of 
interest would not exist in the limit of vanishing momentum transfer because of a massless 
U(l) goldstone contribution. Moreover, while there may exist higher order corrections to 
Eq. ([34]) , this result is exact in the chiral limit at least to the one-loop accuracy. Its 
validity was checked by KZ in the case of supersymmetric QCD, and may also be tested 
in solvable models. 

As the negative sign in Eq. (|33|) is so important for what follows, it is instructive to 
have an alternative derivation of this result, where the sign would be explained in a simple 
and intuitive way. We now present such a derivation. 

The physical meaning of Eq.(^) is very simple: it says that in the chiral limit a 
considerable part of the nucleon mass TEA M p comes from gluons, while the term related 
to quarks (more precisely, to their regulator holds) is equal to — 7 ^- M p . Let us introduce 
chiral combinations of the regulator fields in the standard way: 


3 b 


Rl — 2(1 + 75 )R, Rr — - (1 — 75 )R- 


(35) 


Transformation properties of the R r fields under chiral rotations ~ exp (fays) have 
the very simple form: 


R r Ri —> exp (ia)R r Ri, RiR r —> exp(— ia)RiR r (36) 

The requirement of reparametrization invariance under chiral rotations (^) in the 
limit rn q = 0 can be expressed by equations 

(p\ - M R R r Rie ia \p) = M^e ia p rPl , <p| - = M^e~ ia p lPr , (37) 
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which should be valid for arbitrary a. In this formula AfW is a part of the nucleon mass 
which comes from the fermion regulator fieldQ: = (-^) M p . One can easily check 

this formula once again using an expansion for the scalar density at M R —> oo: 


(p| - M r RR\p) -> (p|^G 2 + 0(1/A4)|p) = M^pp = A M r pp , (38) 


which is exactly the contribution of one flavor to complete expression (^6j). A similar 
calculation for the pseudoscalar part of Eq.fl37l) gives the result (for convenience, we 
multiply both parts of the equation by the factor i): 


(p\ - M R Ri^ 5 R\p) -> (p\^GG + 0(1/M%)\p) = M^pi^p = -^M p pi^p, (39) 


where we used the standard expansion —M R Ri^^R —> |'^GG for the large mass M Rl see 
Eq. (©• Expression ([ID]) is precisely KZ theorem (|TI|). Now we understand the sign in 
KZ formulae ([h|), (]39| ) very well: it has the minus sign because a fermion field gives a 
negative contribution to the nucleon mass in comparison with the main term originating 
from gluons. By the same reasons this term is suppressed in the large N c limit. This 
statement is very clear and unambiguous. 

We would now like to guess a reason why the KZ formula is mostly ignored in the 
literature. The reason is the negative sign of matrix element (]34|), which is absolutely 
crucial in what follows. Had we forgotten about the second term in Eq. (pK|) , we would 
conclude that result ([b|) contradicts the data as it is of a different sign. Our answer to 
this objection is that it is the second charm-induced term in Eq.([T8|) that makes the whole 
expression positive. As will be shown in the next section, an independent estimate of this 
term, based on results on B —> 7 / decays, confirms this claimQ. 

Finally, we are ready to use the KZ theorem given by Eq.([l4|) to find a charm-induced 
contribution to the proton spin from experimental constraint (|l8|). We obtain 


1 

2 M p pi'f 5 p 


(P I 


1 

167r 2 m^ 


g 3 GGG\p) ~ 0.3 


(40) 


We thus conclude that experimental result (|18D , taken together with KZ formula 
requires a substantial contribution of intrinsic charm to the proton spin, which is about 
twice larger than the total singlet light flavor contribution given by the first term in 
Eq. (|THp • We should note that, taken separately, the u and d quark contributions ( |T7| ) 
are numerically considerably larger than the c-quark term. It is the singlet combination 
which is smaller than the charmed quark contribution. Such a smallness of the singlet 
combination seems natural within the large N c approach: it should vanish according 


5 Different n/ fermions multiply this result by the factor n/. 

6 Here we would like to mention that in alternative approach to the anomaly matrix element over 
the proton, which is based on U(l) Goldberger-Treiman relations within a “two-component” approach of 
Ref. Jl|] (see also references in ||), a sign of matrix element (|34|) comes as a result of interplay between 
the r! and a composite gluon ghost field interactions with the nucleon, which may in principle lead to 
any signature. Yet, in the model-independent KZ method this sign is fixed to be negative. A discussion 
of possible consequences of this observation for the formalism of Ref.|l^) would lead us far beyond the 
scope of this paper. 
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to KZ formula (|5ip. We consider this as a reasonable qualitative explanation of our 
result: the charm contribution becomes competitive in comparison with the SU( 3) singlet 
light quark term because the latter is simply suppressed by the factor l/iV c on a natural 
scale. Nevertheless, we understand that this result may be not easy to reconcile with 
one’s intuition. On the other hand, the first term in Eq.(|T8|) is the matrix element of 
a total derivative, i.e. a very unusual object. The fact that it does not vanish can be 
understood as resulting from interactions in a gauge-variant sector of the theory. The 
physical intuition can hardly work in such untypical situation. Thus, to support the 
above conclusion on a role of intrinsic charm in the proton spin, we need to find matrix 
element (f40|) in an independent way. If number (Hi was established without a reference to 
the polarized D1S data, the latter would be completely explained in terms of theoretically 
calculated matrix elements (|T|) and (pUf ) (see Eq.([TS|) above). We will now argue that 
estimate (40) can be conhrmed within a different line or reasoning, based on a very 
different experiment, which is completely independent of both constraint (|T8| ) and KZ 
result fl34p . 


4 Intrinsic charm in the proton and B —> rj data 


We now proceed to an independent estimate of matrix element ([40]) which measures a 
contribution of the intrinsic charm component of the proton spin to the Erst moment 
fdxg\(x) of the polarized structure function gi(x). Somewhat unexpectedly, such an 
estimate can be obtained from new data of B-physics. As will be argued in this section, 
experimental information on B —> rf decays can be used to make conclusions on a role 
of the intrinsic charm component (p|c 7 M 75 c|p) in the proton spin. We will show that this 
quantity may be related to the intrinsic charm component (0| 07 ^ 750 !?/) of the rj. The 
latter matrix element is actually probed experimentally as will be argued below. Thus, in 
what follows we will first spend some time discussing the B-physics, and then demonstrate 
how this can be helpful in the problem of interest. 

Recently, CLEO collaboration has reported 0 results of measurements of inclusive 
and exclusive production of the 7 / in B-decays : 


Br(B -> 7 / + A ; 2.2 GeV < E q , < 2.7 GeV) = (7.5 ± 1.5 ± 1.1) • 10“ 4 , (41) 

Br(B -»■ 7 / + K) = (7.8 ± 1.0) • 10“ 5 . (42) 

Here the inclusive branching ratio contains the acceptance cut intended to reduce a back¬ 
ground from events with charmed meson interactions in a final state. At first sight, the 
above numbers might seem quite innocent. However, simple calculations || [9) reveal 
that these data are in severe contradiction with a standard view of the process at the 
quark level as a decay of the 5-quark into light quarks, which could be naively suggested 
keeping in mind the standard picture of 7 / as a SU(3) singlet meson made of the u —, 
d— and s— quarks. In this picture the B —> 7 / amplitude must be proportional to the 
Cabbibo suppression factor V u b and, as a result, the standard mechanism of the B —> 7 / 
transition yields numbers which are by two orders of magnitude (!) smaller than the data 
for both the inclusive and exclusive cases. Thus, there must be something beyond this 
standard picture. Furthermore, inclusive decays are usually dominated by few-particle 
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final states. Therefore, any meaningful mechanism, intended to explain abnormally large 
numbers (pO] , [42|) , should deal with both the inclusive and exclusive data at the same time. 

A particular mechanism, suggesting a unified description of both the inclusive and 
exclusive modes, was proposed in our recent papers | 8 ], 0. We have shown that at the 
quark level B —> rf decays can be described by the Cabbibo favored b —> ccs process 
followed by a transition of cc into the 7/ . The latter transition is possible due to an 
intrinsic charm component of the r/', and quantitatively can be expressed through the 
matrix element 

(0|c7 M 7 5 c|r/(g)) = . (43) 


A reason why this matrix element may be non-zero can be explained as follows. As the 
c-quark is heavier than the 7 ', it may only exist in the 7/ in a loop. Such a loop with the 
heavy c-quark can be evaluated in terms of gluon fields, populating the ?/, by using the 
background field technique, see Eq. ( 0 )- Therefore, matrix element (|43]) may not vanish 
due to cc gluons transitions. 

Prior to a theoretical calculation of the charm current residue f ( r p into the 7/, it is 

(c) 

instructive to find this quantity “experimentally”. By this we mean a value of /)), needed 
for the above mechanism to explain the data. Simple calculations show that the number 


1 


fjf ~ 140 MeV (“exp”) 


(44) 


yields a good agreement with the data in both the inclusive and exclusive modes. On the 
other hand, this value may look uncomfortably large as it is only a few times smaller than 
the analogously normalized residue ( 0 |c 7 M 7 5 c| 77 c (g)) = if^q^ with f Vc ~ 400 MeV known 
experimentally from the ij c —> 77 decay. Intuitively, these numbers could be expected to 
differ much more drastically as, in contrast to f Vc , the residue ' is a double suppressed 
amplitude. It is Zweig rule-violating and besides contains the 1/m^ suppression factor, 
as it comes from loop effects. Therefore, one could expect it to be very small. In reality 
it is not small. There are two reasons for this. First, m c is not very large on the hadronic 
1 GeV scale. Second, and more important, the Zweig rule itself is badly broken down in 
vacuum CR channels. Of course, such a breakdown contradicts a naive large N c counting 
where a non-diagonal transitions should be suppressed in comparison with diagonal ones. 
However, a more careful analysis |14|, | 8 |] reveals that the large N c picture and breakdown of 
the Zweig rule in fact peacefully co-exist: while the large N c description is quite accurate 
for the 7/, an extent to which the Zweig rule is violated in 7/ yields a large residue We 
stress that the phenomenon of the breakdown of the Zweig rule in vacuum CP channels is 
well known and understood [|T4|], and many phenomenological examples of corresponding 
physics have been discussed in the literature, see e.g. B 0 - The residue (which is 
fundamentally important for our estimates) is another manifestation of the same physics. 

We now proceed to a theoretical analysis of the residue defined by Eq.(f43|). Unfor¬ 
tunately, a detailed theoretical consideration of this quantity would require a repetition 
of our original paper || almost at full length. We thus refer to Ref. || for details, while 
here we would like to give an idea and flavor of our method. It is convenient to start with 
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( 45 ) 


taking the derivative in Eq.TH). Then we obtain 


4' C) = — ( 0 | 2 m c d 7 5 c + . 


' v ' m 2 , 


Since the c-quark is heavy, one can use the Operator Product Expansion in inverse powers 
of the c —quark mass (the heavy quark expansion) 


Of c 


2m c «7 5 c = -—G^G^ - 

47T 


1 


167r 2 m 2 


_,d rabc/^ia s~ib s~ic , 

9 J + ■ ■ ■ 


(46) 


(see the appendix in || for a detailed derivation of this result). Further terms in expansion 
(PH) are neglected in what follows (see however below). We have thus reduced the problem 
to a calculation of the matrix element of the purely gluonic operator: 


Ac) _ 

Jn' ~ 


167r 2 m 2 , m? c 


(OlsV^G" & m CFL W)■ 


(47) 


It may seem at first sight that Eq.([47]) is of little help since matrix elements of gluon 
operators are difficult to calculate. A situation with the r/ is, however, exceptional, as 
the rj' is strongly coupled to gluons and to some extent can be viewed as a remnant of 
imaginary purely gluonic world. In effect, matrix element ( ]-T7| ) is amenable to a theoretical 
study which essentially reduces to a nonperturbative analysis of pure Yang-Mills theory. 
Closely following old ideas due to Witten |T(jj] and Veneziano[|T7J and using some additional 
(nonperturbative) arguments, we have managed to estimate matrix element (|47|) [||: 


A c ) ~ 

J n’ — 


(g 3 G 3 ) 


YM 


47t 2 6 m 2 


Ij^G^G^It]') 


(48) 


Therefore, we have related the residue of the charmed axial current into the 7 / with 
apparently completely unrelated quantity which a cubic gluon condensate in the pure 
Yang-Mills theory (we notice that the matrix element of topological density which appears 
in ([48]) is known {Q\{a s /A^GG^') ~ 0.04 GeV 3 |L4]]). Using all currently available 
information regarding the vacuum condensate ( g 3 G 3 )y-M hr gluodynamics, we have arrived 
at the numerical estimate [HI 


fjf = (50 - 180) Me V 


(49) 


In spite of a large uncertainty of this result, its main source is nevertheless well localized 
and related to a poor knowledge of the cubic condensate in pure gluodynamics. One can 
see that our theoretical estimate ([4U| ) agrees within errors with “experimental” value (fH|) . 

By reasons mentioned above, it is extremely important to calculate matrix element 
(fI7|) in an independent way within a QCD-based model in order to test the general idea. 
Such a calculation based on the instanton-liquid model (see e.g. [|1^] for a review and refer¬ 
ences to original papers) has been carried out |T9| with the following promising result: the 
matrix element of the three-gluon operator as defined by Eq. Q47[ ) is very close to “experi¬ 
mental” value Q44|). We should remind the reader that this model is extremely successful 
in description of low-energy hadronic properties, and a typical accuracy of this approach 
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is not worse than 30 %. A qualitative, microscopic explanation of the enhancement of 


matrix element ([47]) is the following. The model suggests the existence of relatively small 
strongly interacting instantons with a strong gluon field inside the instanton. This field is 
so strong that a relevant parameter gG/m 2 is not small numerically as one could naively 
expect, but rather is of order one. We are quite confident in calculations [[HJ] because a 
similar calculation with the standard gluon current GG is in a good agreement with the 
existing phenomenological data on the rf |L8] . 

We interpret all these results as the evidence that the suggested mechanism indeed 
explains the data on B —* rf given by numbers ((4l| ,|42|) . Therefore, in what follows we 


consider the intrinsic charm component of the rj', described by Eqs. (|43| , [44| , [18| , [49D , as an 
established property of the rf which is probed in experiments and, on the other hand, is 
understood theoretically. Still, in view of a poor accuracy of our theoretical prediction 
(|49|), we will use “experimental” number (|44|) in our subsequent estimates. 

We are now in a position to make our second estimate of an intrinsic charm component 
of the proton spin. Consider the matrix element (q — p' — p, other notations here are 
self-explanatory) 


(N(j/)\cry lt 'y 5 c\N(p)) = u(p') [r^ 5 hi(q 2 ) + q^Mq 2 )] u(p) . 

Taking the derivative and using Eq.(^Bj), we obtain 
1 


(N(p')\ - 


16n 2 m 2 


pGGG\N(p)) = iu(p)q 2 u(p) 2M N h 1 (q 2 ) + q 2 h 2 (q 2 ) 


(50) 


(51) 


Assuming now a rf dominance in this matrix element^] and using the absence of pole 
singularities in the second form factor Ii 2 , which ensures a resolution of the U(l) problem, 
we arrive at the Goldberger-Trciman type relation 


MO) — „ n r SrfNNflf 


2 M; 


(52) 


N 


or, equivalently, 






(53) 


Note that owing to Eq.(j46j) the anomaly term cancels out in (|53|)fj. Unfortunately, the 
precise value of (] v >nn is not known, and phenomenological estimates of the coupling g v 'NN 

7 Note that such a saturation becomes exact in the large N c limit. 

8 This is in contrast to the standard case with U(l) Goldberger-Treiman relations for the SU (3) singlet 
light flavor current discussed in Ref. 0- The anomaly term is a total derivative, and the fact that it gives 
a non-zero contribution to the proton matrix element is due to the presence of non-pole subtraction terms 
in dispersion relations, which have the same nature as subtraction terms in a correlation function of the 
topological density 0 , and can be assigned to a ghost contribution. This is the reason why the proton 
spin has a “two-component” form in the approach of Ref.[0. As a result, the light quark component of 
the proton spin is expressed in terms of two unknown quantities, neither of which is directly measurable. 
For these reasons we prefer to use the exact KZ theorem for the anomaly matrix element rather than 
results from 0. On the contrary to the case [0, the anomaly term does not show up in Eq.(0), and 
this is why it is the physical coupling constant g v 'NN that stands there. 
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vary from g^NN — 3 to g^NN — 7 [PU| . Assuming the same range for this quantity in the 
chiral limit and using numerical value (fUD, we arrive at the estimate 


1 

2M n u{jp)i^u{p) 


{N(p) I 


T J^g*GGG\N(p)) = 0.2 - 0.5 . 


(54) 


This result has a large uncertainty mainly due to a poor knowledge of g v 'NN- In view of 
a large uncertainty in g v 'NN, we have neglected a perturbative evolution of 

A few remarks are in order. First, as we already mentioned, the microscopic pic¬ 
ture suggests that the matrix element under consideration is not small because of gluon 
fluctuations which could be large. In different words, the parameter of the heavy quark 
expansion gG/m 2 c could be of order one. Such an explanation might imply a bad conver¬ 
gence of expansion (|46|) where we keep the first term and neglect all the rest. A scientific 


answer to the question of convergence of series (|46l) would require the knowledge of higher 
dimensional matrix elements which are not available at the moment. Alternatively, one 
could try to address this issue within the instanton liquid approach [jnj. ffowever, irre¬ 
spectively to an outcome of such an analysis, we could consider “experimental” value of 
(44), which effectively takes into account all powers of the 1 jm c expansion. Therefore, 


all of them are also effectively included in the left hand side of Eq. (^4[) . Thus, a possible 
bad convergence of expansion (f|6|) does not affect our results at all. ffence we see that 
our final estimate ( HI ) confirms the above conclusion on a large c-quark component in the 
proton spin, as its sign and order of magnitude agree with phenomenological constraint 
(|40j). The large magnitude of the charmed current residue f into the rj' is the main fac¬ 
tor giving rise to large number (|5|). We emphasize again that estimate (|54|) is obtained 
in an entirely independent of the DIS phenomenology and KZ theorem way, and relies on 
the data and theoretical results of B-physics. If we fix the r/NN coupling constant in the 
chiral limit by “experimental” constraint (|4Q|) , it is consistent within uncertainties with 
independent theoretical prediction (|5|). We thus conclude that the polarized DIS data 
agree within errors with theoretical results expressed by Eqs.( 


D, (fp and (HI). This 
implies that the main contribution to T)' is due to the charm component of the proton 
spin. We now proceed to a discussion of possible ways to test this conclusion. 


5 Signals for intrinsic charm in polarized experiments 

In this section we would like to argue that our explanation of the polarized DIS data 
has very dehnite consequences for future polarized experiments which thus will be able 
to test the picture suggested by this paper. We will only deal with qualitative aspects 
of expected phenomena, which follow immediately from our results. A more detailed 
analysis is beyond the scope of this paper and deserves a separate study. 

We have found that the light quark and charm contributions to the first moment Tf 
have different signs and magnitudes. Here and in what follows any references to the proton 
spin are understood in the sense of the singlet part of the spin operator AS + 2Ac only. 
While the former is negative, the latter is positive with about a twice larger magnitude. 
The immediate consequence of this fact is that an integral asymmetry for events with no 
charm in the final state is expected to be negative with approximately the same magnitude 
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as a total integral asymmetry. This prediction may hopefully be tested experimentally, 
provided the two types of final states can be distinguished. 

Yet, a most interesting test of our scenario is related to a charm distribution in charmed 
final states. We find that the charm contribution to Tf is about an order of magnitude 
larger than results obtained within perturbative schemes |^T] . The difference from previ¬ 


ous analyses comes due to strong nonperturbative fluctuations in the axial SU(3) singlet 
channel, which lead to a large magnitude of the intrinsic charm component of the pro¬ 
ton spin. Our definition of the latter quantity matches the well known definition |^2| of 
the intrinsic charm in terms of a Fock state description as a multiple connected loop of 
the c-quark in the protonj]. In our approach this property is used in a constructive way. 
Furthermore, one can expect that the intrinsic charm component of the proton spin will 
be traced in polarized experiments. We note that in the case of unpolarized DIS various 
mechanisms, which are able to identify the intrinsic charm in the proton, have been dis¬ 
cussed in the literature. In particular, they include a liberation of the charm as a result of 
scattering on light quarks ]^4|, |25|, or direct DIS on the intrinsic charm JTT]] . The intrinsic 


charm contribution to the proton spin, considered in this paper, corresponds to the latter 
case. 

There exists a clear distinction between pictures of a charm distribution in the final 
state corresponding to different mechanisms of the charm production in DIS. A vast 
literature |26| is devoted to the charm production in polarized experiments as an effective 
tool to measure the gluon spin dependent distribution A G(x, Q 2 ). This proposal is based 
on the perturbative photon-gluon fusion (PGF) mechanism which assigns a production 
of the cc pair to the hard subprocess 7 *g —> cc. It is undoubtably true that the PGF 
process adequately describes a heavy quark production in the case of usual unpolarized 
DIS, and expresses it in terms of the unpolarized gluon distribution. Moreover, in the 
unpolarized case the manifestation of the intrinsic charm component of the nucleon is 
still an open problem (ll|, ^5[. We believe that the underlying reason for a relatively 
subdominant role jT7| of the intrinsic charm in usual DIS is a strong suppression of 
nonperturbative effects in vector channels. However, as we argued in Introduction, it 
would be potentially dangerous to automatically transfer concepts and results, valid for 
unpolarized DIS, to the polarized case. Although at the perturbative level results look 
similar, these arguments ignore important differences between the two cases, which only 
show up beyond perturbation theory. That this way of thinking is fallacious can be 
most clearly illustrated on the example of the pseudoscalar nonet. The Zweig rule is 100 
% violated there by nonperturbative effects, though no difference from the Zweig rule 
conserving case of vector mesons is seen in perturbation theory. 

As, in contrast to the perturbative PGF mechanism, polarized DIS on the intrinsic 
charm does not have a strong a s suppression, we expect that the charm production in 
polarized experiments will be overwhelmed by events due to intrinsic charm mechanisms. 
Quantitative predictions in this case are troublesome as, in addition to direct DIS on 
the intrinsic charm, there exist other possibilities for the extrinsic and intrinsic charm 
production in DIS [F~ 


25]. Model independent methods needed for their evaluation 


are not available at present, though estimates with model light cone wave functions for 


3 See |2li|] for a recent discussion of the intrinsic charm physics in a different content. 
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the intrinsic charm [^] are possible. Nevertheless, different mechanisms for the charm 
production can be disentangled on the basis of a event topology analysis |TT], The 
PGF mechanism corresponds to two jets (plus a spectator) events with a charmed hadron 
(typically a D-meson) carrying a large fraction of the photon energy in the current jet 
and an anticharm (e.g., D ) hadron in the target jet. A topology of these PGF events is 
expected to deviate from the planar one [|ll|]. The cc pair has a low invariant mass as 
the fractional energy of the proton carried by the gluon is typically small. The transverse 
momentum and invariant mass of the pair will grow with the photon virtuality Q 2 . As 
has been argued above, we expect that these events will be parametrically suppressed by 
powers of a s in comparison to events reviving the intrinsic charm in the nucleon. 

As for DIS on the intrinsic charm, one can distinguish between two types of processes. 
The first one is an indirect process when the photon scatters on a light quark with a 
subsequent liberation of the charm 


24, 


. In this case both the charm quark and 
antiquark are in the target fragmentation region, and can hadronize into both open charm 
(e.g. D,A C ) or hidden charm {J/ip, etc.) hadrons. In this process, the intrinsic charm 
shows up through a nonperturbative final state interaction. Charmonium states, produced 
in this region, will presumably have a large fraction of the proton momentum with pj_ ~ m c 
and independent of Q 2 . On general grounds, one may expect an excessive production of 
S-wave quarkonia (e.g. r ] c ) in comparison to P-wave states. 

The second process is direct DIS on the intrinsic charm quark [III]. The target charm 


quark is emitted at large forward rapidities, while the scattered quark has somewhat 
smaller near forward rapidities. The difference in rapidities grows with Q 2 . The average 
transverse momentum of the cc pair is not expected to grow with Q 2 . Thus, our results 
allow one to suggest that a substantial (presumably dominant) part of charm events in 
polarized experiments at large Q 2 will be an open charm hadrons (or S-wave quarkonia) 
production with low p± ~ m c . 

Finally, we note that all said above on the charm production can be extended to beauty 
production. The latter is suppressed relatively to the former by the factor ml/m 2 ~ 10, 
which may imply that the beauty can be produced nonperturbatively via the intrinsic 
beauty mechanism at approximately the same level as the charm production through the 
perturbative PGF process. 


6 Conclusions 

In this paper we have presented a somewhat unorthodox point of view of the data on 
the first moment fdxg±(x ) of the spin g±(x) structure function measured in polarized 
DIS experiments. We have abandoned attempts to explain the data within the quark 
model or related to it held theoretical models, and reformulated the problem in terms 
of particular matrix elements of quark and gluon operators. In our opinion, none of the 
phenomena, which are of crucial importance for understanding the proton spin problem 
(spontaneous breaking of the chiral invariance, a resolution of the U(l) problem, the 
dimensional transmutation, etc.), can be adequately described within the parton or quark 
model. On the contrary, all of them are essentially nonperturbative, and need to be 
addressed within nonperturbative QCD. We have further shown that the data can be 
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interpreted as a constraint on the sum of matrix elements of the anomaly and intrinsic 
charm operators. The latter object is a new principal element of our analysis, which has 
not been considered hitherto in the literature in the content of the proton spin problem. 
Still, we have argued that the intrinsic charm component of the proton spin is quite large 
and constitutes a main contribution to the first moment Tf. Its large magnitude is related 
to strong nonperturbative fluctuations in the axial SU (3) singlet channel. Moreover, the 
inclusion of the intrinsic charm proton spin in the analysis is absolutely necessary as it 
is the only way to reconcile the polarized DIS data with both an exact Kuhn-Zakharov 
low energy theorem on one side, and new data and theoretical results on B —> rj f decays 
on the other. We believe that these seemingly uncorrelated problems are actually tightly 
connected. The situation reminds us the J/-0 discovery in 1974, when a charmonium 
state (“hidden charm”) was observed simultaneously in e + e _ collisions at SLAC [|25| and 
at the proton machine at Brookhaven |29|. We believe that we are now facing a similar 
case, when different experimental groups see the “intrinsic charm” in polarized DIS ( ]T0| ) 
and B-decays (HO simultaneously. 

We feel that at a microscopic level, the new effect of large intrinsic charm fluctuations 
originates in the instanton physics [|18|], where strong gluon fields in the singlet channel 
are able to lead to very unexpected phenomena. Our results seem to imply a kind of 
universality of this nonperturbative physics with universality classes defined by quantum 
numbers relevant for a process considered, but not concrete particles involved. 

Our explanation of the polarized DIS data has very definite consequences for the future 
polarized experiments, and thus can be tested there. 
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